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This Letter studies quantum limits to dynamical sensors in the presence ol decoherence. A 
modified purification approach is used to obtain tighter quantum detection and estimation error 
bounds for optomechanical force sensing and optical phase sensing. When optical loss is present, 
these bounds are found to obey shot-noise scalings for arbitrary quantum states of light under 
certain realistic conditions, thus ruling out the possibility of asymptotic Heisenberg error scalings 
with respect to the average photon flux under those conditions. The proposed bounds are expected 
to be approachable using current quantum optics technology. 
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The laws of quantum mechanics impose fundamental lim- 
itations to the accuracy of measurements, and a fun- 
damental question in quantum measurement theory is 
how such limitations affect precision sensing applications, 
such as gravitational-wave detection, optical interferom- 
etry, and atomic magnetometry and gyroscopy [1] . With 
the rapid recent advance in quantum optomechanics Q 
and atomic 3] technologies, quantum sensing limits have 
received renewed interest and are expected to play a key 
role in future precision measurement applications. 

Many realistic sensors, such as gravitational-wave 
detectors, perform continuous measurements of time- 
varying signals (commonly called waveforms). For such 
sensors, a quantum Cramer-Rao bound (QCRB) for 
waveform estimation (3 and a quantum fidelity bound for 
waveform detection [5[ have recently been proved, gen- 
eralizing earlier seminal results by Helstrom Q. These 
bounds are not expected to be tight when decoherence 
is significant, however, as Refs. use a purification 

approach that does not account for the inaccessibility 
of the environment. Given the ubiquity of decoherence 
in quantum experiments, the relevance of the bounds to 
practical situations may be questioned. 

One way to account for decoherence is to employ the 
concepts of mixed states, effects, and operations [7(- Such 
an approach has been successful in the study of single- 
parameter estimation problems 0, Q, but becomes in- 
tractable for nontrivial quantum dynamics. To retain 
the convenience of a larger Hilbert space, here I extend 
a modified purification approach due to Escher et al. 
@ and apply it to more general open-system detection 
and estimation problems beyond the paradigm of single- 
parameter estimation considered by previous work [1, Q . 
In particular, I show that 

1. A quantum model of optomechanical force sens- 
ing can be transformed to an optical phase sens- 
ing problem with classical phase shift, such that a 
unified formalism can treat both problems andpro- 
duce tighter bounds than the results in Refs. 0, 



2. For optical phase waveform detection with loss and 
vacuum noise, the errors obey lower bounds that 
scale with the average photon flux akin to reduced 
shot-noise limits, provided that the phase shift or 
the quantum efficiency is small enough (the pre- 
cise conditions will be given later). This rules out 
Heisenberg scaling of the detectable phase shift 
[loj in the high-flux limit under such conditions, 
as well as any significant enhancement of the er- 
ror exponent by quantum illumination jllfl in the 
low-efficiency limit. 

3. The mean-square error for lossy optical phase 
waveform estimation also observes a limit with 
shot-noise scaling, which generalizes the single- 
parameter results in Refs. y|, Q and rules out the 
kind of quantum-enhanced scalings suggested by 
Refs. 12J in the high-flux limit. 



These results not only provide more general and realis- 
tic quantum limits that can be approached using current 
quantum optics technology [l3, 14 1, but may also be rel- 
evant to more general studies of quantum metrology and 
quantum information, such as quantum speed limits [IH 
and quantum phase transitions 

Let x be the vector of unknown parameters to be es- 
timated and y be the observation. Within the purifi- 
cation approach 0, HJ , the dynamics of a quantum sen- 
sor is modeled by unitary evolution (U x as a function 
of x) of an initial pure density operator p = 
and measurements are modeled by a final-time positive 
operator-valued measure (POVM) E(y) using the princi- 
ple of deferred measurement [17[ . The likelihood function 
becomes P(y\x) = tr [E(y)U x pU^\ . For continous-time 
problems, discrete time is first assumed and the continu- 
ous limit is taken at the end of calculations. Refs. [1, [B| 
derive quantum bounds by considering the density oper- 
ator U x pU\. 

Suppose that the Hilbert space {H = Ha ® "Hb) is 
divided into an accessible part (Ha) and an inaccessi- 
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FIG. 1. Quantum circuit diagram [17J for the modified purifi- 
cation approach. 



ble part (Hb)- The POVM should now be written as 
E(y) = E A {y)® l B , where E A (y) is a POVM on H A and 
1 b is the identity operator on Hb, which accounts for 
the fact that Hb cannot be measured. The key to the 
modified purification approach, as illustrated by Fig. [TJ is 
to recognize that the likelihood function is unchanged if 
any arbitrary x-dependent unitary Ub on Hb is applied 
before the POVM: 



where 



P(y\x)=tr{[E A (y)®l B }p x {U B )}, (1) 



p x (U B ) = (l A ® Ul)U x pUl{l A ® 17 B ) (2) 



is a purification of tiB(U x pU^.), such that ItbPx = 
^B{U x pU\). Judicious choices of Ub can result in tighter 
quantum bounds as a function of p x (Us)- 

First, suppose that x^ and x^ are two hypotheses for 
x and x(y) is the estimate. The following theorems are 
applications of the modified purification and Helstrom's 
bounds for pure states @: 

Theorem 1 (Fidelity bound, Neyman-Pearson crite- 
rion). For any POVM measurement E A (y) of tvs p x 
in H A , the miss probability Pqi = J x=0 dyP(y\x = 1) 
given a constraint on the false-alarm probability Pig = 
J-—i dyP(y\x = 0) < a satisfies 



Poi > P(a,F) 

= f 1- [Vc7F + y/(l-a)(l 
~ I 0, 



F)Y 



a < F, 
a> F, 



(3) 



where F is the fidelity between the following pure states 
in H A ®Hb: 



po = UoWWUl p x 
F(po,Pi) 



U B Ui\^)^\u\u B , (4) 
V\ululU \V) \ (5) 



p m = p x ( m ) , U m = U x {m), and 1 a ®Ub, abbreviated as 
Ub , is an arbitrary unitary onHs- 

Theorem 2 (Fidelity bound, Bayes criterion). The av- 
erage error probability P e = P\qPo + Pq\Pi with prior 
probabilities Pq and Pi = 1 — Pq satisfies 



Pe > 



i(W 



- Jl-±P PxF \ > PaPxF. 



(6) 



Proof of Theorem[]]and{^ Ref. p shows that the 
bounds with the likelihood function P(y\x) ~ tr[E(y)p x ] 
are valid for any POVM E(y) on H A ® Hb, so they 
must also be valid with P(y\x) — tx A [E A (y) txs p x ] — 
tr{[E A {y) ® l B ]p x } for any POVM E A (y) on H A . □ 

Since the lower bounds are valid for any Ub, Ub 
should be chosen to increase F and tighten the bounds. 
The maximum F becomes the Uhlmann fidelity between 
mixed states tiB po and trs pi 17, 3]. The bound on P e 
obtained using this method is thus weaker than the Hel- 
strom bound for the mixed states [(| , although it can be 
shown that the error exponent for the Uhlmann-fidelity 
bound is within 3dB of the optimal value (l9| . 

Next, consider the estimation of continuous parameters 
x with prior distribution P(x). A lower error bound is 
given by the following: 

Theorem 3 (Bayesian quantum Cramer- Rao bound). 
The error covariance matrix S = E (x — x) (x — x) T = 
J dydxP(y\x)P(x) (x — x) (x — x) satisfies a matrix in- 
equality given by 



where 



J 



(Q) 



j(C) 
> 3 k 



4 / dxP{x) 



d 2 F(p x <,p x ) 



dx'ftx' k 



dxP(x) 



dlnP(x) d\nP{x) 



dxj 



dx k 



(7) 

(8) 
(9) 



Proof. See Ref. [J] for a proof and Refs. [2(J] for the rela- 
tion between J^> and F. □ 

Ub should again be chosen to reduce the quantum 
Fisher information (QFI) J^\Ub) and thus tighten the 
QCRB. 

As an application of the general theory, consider the 
estimation of a force x(t), t € [to, tf], on a quantum mov- 
ing mirror via continuous optical measurements. For sim- 
plicity, assume that any optical cavity dynamics can be 
adiabatically eliminated. Let (q,p) be the mechanical 
position and momentum operators, and a(t) be the an- 
nihilation operator for the one-dimensional optical field, 
which obeys [a(t), a^(t')] = S(t - t') Suppose 



U x = Texp 



1 r*f 

- / dtH(t) 



with a Hamiltonian given by 

H = H B (q,p,x(t)) - 2hMkqI(t), 



(10) 



(11) 



where T is the time-ordering superoperator, Hb is the 
mechanical Hamiltonian, M is the effective number of op- 
tical reflections by the mirror, k is the optical wavenum- 
ber, and I(t) = (t)a(t) is the photon flux. 
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In an optomechanics experiment, the mechanical oscil- 
lator is measured only through the optical field, so one 
can let be the mechanical Hilbert space. Assume 



u B(tf) = exp 



dtH B (q,p,x(t)) 



to 



(12) 



and let Uj = U B (tf)U x , which can be calculated using 
the interaction picture. The result is 



Ui = Texp 



2iMk I dtqi(t)I(t) 

to 



(13) 



with qi(t) = U B {t)qU B {t) and p x = UipU\. If the 
mechanical dynamics is linear, qi(t) can be expressed 
in terms of the mechanical impulse-response function 
h{t,t') as q r (t) = q (q,p,t) + dt'h{t,t')x{t'). Apart 
from the transient solution qo, the mechanical Hilbert 
space has been removed from the model using U B , and 
the problem has been transformed to a problem of op- 
tical phase sensing, with the phase given by <f>(t) — 
2Mk Jl f dt'h(t,t')x(t'). 




FIG. 2. A model of the lossy optical phase sensing problem. 

I now turn to the problem of optical phase sensing, as 
depicted in Fig. [2] Consider first the detection problem. 
Let cj>(t) be the phase shift between the two hypotheses, 
Uq = UaUab, and U± — Uab, where 

ftf 

U A = exp i dt(j>(t)I{t) , (14) 

I J to 

Uab = exp jm J f dt [a f (t)b(t) + a(t)b ] '(*)] j. (15) 

Uab models loss as a beam-splitter coupling with an- 
other optical field b{t) in vacuum state \Q) B before the 
phase modulation, such that |\E f ) = \iP)a®\0)b- Uab can 
also account for loss after the modulation, as shown in 
Ref. [22j. The fidelity becomes 



With 



F= (V\tf AB U B U A U AB m 



U B = exp i I dt\(t)<p(t)tf (t)b(t) 



(16) 



(17) 



where A(t) is a free parameter to be specified later, 
Eq. (fl"6)) can be simplified usin g th e SU (2) disentangling 
theorem [2^], as shown in Ref. |22| . The result is 



F 

z{t) 







exp 


[/' 

I J t 



dtl(t)lnz(t) 
) + (l-77)e- iA WW, 



(18) 
(19) 



with r\ = cos 2 k defined as the quantum efficiency. For 
example, the coherent-state value is 



F coh = exp 



477 / dt(I(t)} sin 



2 



(20) 



where (O) = (ip\0\ip). Measurements that can saturate 
the Bayes error bound in Theorem [5] for coherent states 
are known 0, [l3| . 

To bound F in general, Jensen's inequality can be used 
if z is real and positive. The following lemma provides 
the necessary and sufficient condition: 

Lemma 1. There exists a A such that z = rje 1 ^ + (1 — 
7y)e _ ' A ^ is real and positive if and only if one of the fol- 
lowing conditions is satisfied: 



(I) : V < 



1 

2' 



(21) 



1 1 - n 

(II) : j] > - and | sin </>| < and cos <f> > 0. (22) 

2 77 

Proof. Consider the circle traced by z{\) centered at r/e 1 ^ 
with radius 1 — n on the complex plane, z = \z\ > for 
some A is equivalent to the condition that the circle inter- 
sects the positive real axis, for which the necessary and 
sufficient condition is given by one of Eq. (|2"Tj) (the circle 
encloses the origin for any (j> and thus always intersects 
the axis) and Eq. (|2"21 (the circle intersects the axis for 
some <f> on the right-hand plane only). □ 

Eq. (|2"2"|) holds when <f> is sufficiently small. For ex- 
ample, M ~ 100, |g| ~ 10" 19 m, 2?r/fc ~ 1 /im, 
and (1 - rj)/<n - 10~ 2 for LIGO 0, leading to ~ 
2Mk\q\ - 10 -10 , and Eq. {22} is easily satisfied. 

The following theorem is a key technical result of this 
Letter: 

Theorem 4. If Eq. {Hp or is satisfied for all 4>(t), 



F > exp 



dt(I(t)) In \z(t){ 



to 



F z . (23) 



Proof. With z = \z\ > under the condition in 
Lemma Q] and writing as a superposition of eigen- 
states of I(t), one can apply Jensen's inequality and get 
(exp(/ dtlhxz)) > exp(J dt(I)\nz). Eq. ^ then fol- 
lows from Eq. ([Tg]). □ 
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Compared with Eq. (|20l) . F z has the same shot-noise 
scaling with respect to the average photon flux (7), as 
— In F z scales linearly with (I) . Since error-free detection 
with F = is possible with pure states (25|, this shot- 
noise-scaling bound is a very strong result. It should 
also have implications for M-ary phase discrimination in 
general [26j |. 

The following corollaries are some analytic conse- 
quences of Theorem 0] that exemplify its tightness: 



Corollary 1. If |0| < 1 and 1 770/(1 - i])\ < I, 



F z « exp 



V 



VJt 



dt{i{t))4, 2 {t) 



(24) 



which determines the exponent of F z in Eq. (123 
leads to Eq 



Proof. Let Ao = 77/(1 — 77). Since Imz = 77 sin — (1 — 
7 7 )sinA0 = 0, A = A O [1 + O(0 2 ) + O(A 2 2 )], \z\ 2 = 1 - 
47/(1-7?) sin 2 [(l + A)0/2] = 1-A O 2 [1 + O(0 2 ) + O(A 2 2 )], 
and ln|z| 2 = -A O 2 [1 + 0{(j> 2 ) + 0(A O 2 ) + O(\ 2 (j) 2 )}, 

and 
□ 

Eq. (HU) differs from F coh « exp(-77 J dt{I)4> 2 ) by just 
a constant factor of 1/(1 — 77) in the exponent. 

To obtain another measure of detection error, I formal- 
ize the concept of detectable phase shift [l(| as follows: 

Definition 1 (Detectable phase shift). A detectable 
phase shift 4>'(t) given acceptable error probabilities a' 
and f3' is a <p(t) that makes P\q < a' and Pq± < ft' . 

Corollary 2. Under the conditions of Corollary]^ 



dt{i{tM 2 {t) 



> 



i{-hxF') 



(25) 



where F' = vaaxpi>pr a ' t F) F with (3 defined in Eq. f3|). 

Proof. Any achievable (P10, Pqi) must lie above the con- 
vex curve Pqi = (3(Pi ,F). This means that for P 10 < a' 
and Pqi < (3', f3' > /?(«', F) must hold, and hence 
F < max^>^ (a , )F) F = F'{a',[3'). Since F > F z , 
F z < F' must hold for the constraints on (P10, P01) to be 
possible. Eq. ([25|) then follows from Eq. j24|). □ 



For example, if (I(t)} is constant, 

1 C l f I -- /; 



, , , dt^(t)^ 



(- InP'), 



(26) 



with N = (tf — ti)(I). This is lower than the shot-noise 
limit by a constant factor of 1 — 7/ only, rulin g ou t the kind 
of Heisenberg scaling suggested by Refs. [l0( for lossy 
weak phase detection in the ./V — > 00 limit. 



Corollary 3. //i)«l,F»Rif, 



coh ■ 



Proof. Since Imz = 0, = sin 1 [77 sin 0/(1 — 77)] = 
0(77), z = Re z = 77 cos (f> + (1 — 77) cos \cf> = 1 + 77 cos <fi — 
77 + 0(?7 2 ), and ln|z| 2 = -4?7 sm 2 ((f>/2) + 0(rf), which 
leads to - In F z = -[1 + 0(rj)] In F coh . □ 



Corollary [3] proves that, analogous to the case of target 
detection [27] , the coherent state is near-optimal for any 
phase detection problem in the low-efficiency limit, ruling 
out any significant enhancement of the error exponent by 
quantum illumination llj in this scenario. It remains an 
open question whether quantum illumination is useful for 
high-thermal-noise low-efficiency phase detection. 

Consider now the estimation problem. The QFI 
jW(tj,t k ) ^Jhnst^oJ^/St 2 for estimating <p(t) 



found in Ref. |22j to be 



is 



J^(t t') 

* I =[V- A(l - V)} 2 (A/(t)AJ(f )> 

+ (1 + A) 2 77(1 - 77) (/(iMi-i'), (27) 

where A/(/j) = I{t) — (I(t)) and A is assumed to be con- 
stant for simplicity. If the light source has stationary 
statistics, {I it)) is constant, (AI(t)AI(t')) depends on 
t—t' only, and a power spectral density Sai(u) can be de- 
fined by {AI(t)AI(t')) = j^idu /2ir)S AI (uj) exp[icuit - 

if)]. A spectral form of jf ] is then J { ®\uj)/A = 
[7)-A(1-77)] 2 ,SajH + (1 + A) 2 7)(1 -77)(/). The minimum 
QFI with respect to A is hence 



minjf ) H=4 



Sai(w) 77(7) 



(28) 



This method of deriving an open-system QFI is the same 
as the one proposed in Refs. [9j and the result is a gener- 
alization of earlier ones for lossy static-phase estimation 
in Refs. For <j>{t) = f™ x dt' 'g(t-t')x(t'), the QCRB 

on E t ee E[i{t) ~ xit)] 2 becomes 0| 



St > 



00 2tt | fl ( w )|2 m in A J^H + Ji^(w) 



CO), 



AO, 



(29) 



where g(u>) is the Fourier transform of git) and Jx (ui) 
is the prior information in spectral form. Compared 
with the coherent-state value 7^| )h (a;) = 4/7(7) , the QFI 
for any state is limited by the same shot-noise scaling 
(min^ j{ (w) < 47/(7)/(l — 77)). This rules out the kind 
of quantum-enhanced scalings suggested by Refs. [12;] in 
the high- flux limit when loss is present. 

For optomechanical force sensing, the detection and es- 
timation error bounds here are more realistic to achieve 
than the bounds in Refs. 0, IH, as the former should 
be approachable using the quantum optics technology 
demonstrated in Refs. LJ, [lj], whereas the latter also 
requires measurement backaction noise to dominate the 
fluctuation of q and may require quantum noise cancel- 
lation (US]. 

Discussions with R. Nair, H. Wiseman, and C. Caves 
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the Singapore National Research Foundation under NRF 
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Order of loss and phase modulation 



Here I prove that the reduced state tr B p x is the same regardless of the order of the optical loss Uab and the phase 
modulation U A , viz., 

Lemma 2. 

tr B (u A U AB pU AB U A ) = tr B (u AB U AP U A U AB ) (30) 
if p = p A ® Pb and p B is a thermal state. 
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Proof. It suffices to prove it for one mode in H,a and one mode in Us', generalization to the multimode case is 
straightforward. Suppose 

U A = exp(i0a^a), Uab = exp[i/i(a^6 + ab' )], Ub = exp^tpb^b), (31) 

where [a, at] = [6,^] = 1. Then 

[/A^Asf/l = exp[iK(e^a t 6 + e"^a& + )] = U B U AB U B , (32) 
Was = U B U AB U B U A . (33) 

Since UbPbU b = pb for a thermal state, 

tr B (VaC/ab^as^) = tr s (u b UabU b UapU a uIu ab U b ) (34) 
= tr B (u^UabUapU^U^Ub^ (35) 
= tr B (l/Afll/ApIT-t l/j^) . (36) 

□ 

With the concatenation property of thermal-noise channels [29| , any optical loss with thermal noise at any stage of 
a phase modulation experiment can be modeled by a single beam splitter before or after the modulation. 



SU(2) algebra 

Consider again the single-mode case with 

U B = exp(iA06 t 6). (37) 

First, compute the following quantity using the Heisenberg picture: 

U ab UbUaUab = exp(i#), g = a'V - \b%' , (38) 

a' = cos na + i sin Kb, b' = cos nb + isinna. (39) 



This gives 



g = pa) a + vb^b + ij(a^b - atf), (40) 

p = cos 2 k — A sin 2 k, (41) 

v = sin 2 k — A cos 2 K, (42) 

7 = (1 + A) sin k cos K. (43) 



Next, define SU(2) operators as 



J_=a)b, J + = ab\ J 3 = - (tfb - a) a) , J = - (tfb + a 1 a) , (44) 

where J commutes with the rest of the operators. In terms of the redefined operators, 

exp(i<70) = exp [i<fi(p + v) J] exp [i(j)(v — p)Js — 7</>J- + l<t>J+] ■ (45) 
The following theorem is useful: 

Theorem 5 (SU(2) disentangling theorem). Given J± and J3 that obey the commutation relations 

[J+,J_]=2J 3 , [Js,J±] = ±J±, (46) 

exp (A+ J+ + A_ J_ + A 3 J 3 ) = exp (A+ J+) A3 3 exp (A_ J J) , (47) 

where 
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Proof. See, for example, Ref. [22 
For the case of interest here, 



A 3 = i<f>(v-[i), X+=-f<j), A_ = --yep, 

1 + A, 



A, 



i + A , , „ . . i + A ; 

cos — - — <p — i[l — 2?]) sm — - — tp 



The disentangling theorem is useful because exp(A_ J_)|0)s = \0)b- 



| exp(i50)|O) B = B (0|exp [i<fi([i + v)J] exp (A+ J+) A3 3 exp (A_ J_) |0) 
= exp [«<^(/i + u)a ! a/2\ A3 a ^ 2 



z = -qe^ + (1 - i])e 



-iXcf) 



□ 

(49) 
(50) 

(51) 

(52) 
(53) 
(54) 
(55) 



Generalization to the multimode case is straightforward. For continuous optical fields, a(t) can be first discretized in 
time as a(tj) « \/Siaj before applying the multimode result and taking the continuous limit. 



Quantum Fisher information 

Consider the multimode case with annihilation operators Oj and bj. Let 



U A = exp I i^foa] 



Uab = exp 



Ub = exp f iX <t>jbjb. 



The QFI matrix jW) can be computed by considering the fidelity for small <j>j and F w 1 [2C 



F = 



*|C^ B 4C/ A C/ AB |*) 2 « 1- i^^.jW)0 fc + O(||0|| 4 



(56) 



(57) 



This also shows why F is more difficult to calculate than in general, as j'^J is just a second-order term in F. 
Writing the fidelity as 



*|exp N^ft^-W 



I get 



After some algebra, 



J 



3j = fia^dj + + i"f(ajbj — ctjbj), 
- 4(*|A 5 ,A.9fc A//, = ft - (*| fe |*>. 



t(O) 

7 ifc 



= (cos 2 k - A sin 2 k) 2 (AnjAn k ) + (1 + A) 2 sin 2 n cos 2 k (nA S jk , 



(58) 
(59) 

(60) 
(61) 



where rij = ajaj. In the continuous-time limit with tj = < + 



aj — > a(tj-), 



— nj -> I(tj), 



j(Q) 

-^^j(Q)( tj ,t k ), (62) 



and Eq. (27) in the main text is obtained. 



